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We investigate growing networks based on Barabasi and Albert's algorithm for generating scale- 
free networks, but with edges sensitive to overload breakdown. The load is defined through edge 
betweenness centrality. We focus on the situation where the average number of connections per 
vertex is, as the number of vertices, linearly increasing in time. After an initial stage of growth, 
the network undergoes avalanching breakdowns to a fragmented state from which it never recovers. 
This breakdown is much less violent if the growth is by random rather than preferential attachment 
(as defines the Barabasi and Albert model). We briefiy discuss the case where the average number 
of connections per vertex is constant. In this case no breakdown avalanches occur. Implications to 
the growth of real-world communication networks are discussed. 
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I. INTRODUCTION 



II. DEFINITIONS 



Large sparse networks are the underlying structure 
for transportation or communication systems, both man 
made (like computer networks or power grids |^) 
or natural (like neural networks |^ or biochemical net- 
works iQ). These networks displays both randomness 
and some self-induced structure influencing the flow of 
transport and robustness against congestion or break- 
down in the network. One of the most conspicuous 
structures among real-world communication networks is 
a highly skewed distribution of the degree (the number 
of neighbors of a vertex) [0, ^ ^. 



Avalanching breakdown in networks where the edges 
or vertices are sensitive to overload is a serious threat 
for real- world networks. A recent example being the 
black-out of 11 US states and two Canadian provinces 
the 10th August 1996 0. Recently the overload break- 
down problem for vertices in growing networks with an 
emerging power-law degree distribution has been stud- 
ied S. In the present paper we investigate the overload 
breakdown problem when edges (rather than vertices) are 
sensitive to overloading. We use the standard model for 
such networks — the Barabasi- Albert (BA) model p^ , 
but with a maximum load capacity assigned to each 
edge. The load is defined by means of the betweenness 
centrality — a centrality measure for communication and 
transport flow in a network ||ll[| . The procedure enables 
us to study overload breakdown triggered by the redistri- 
bution (and increase) of load in a growing network. This 
is in contrast to earlier models of cascading breakdown 
phenomena, all dealing with vertex breakdown, that has 
taken a fixed network as their starting point |12, 0|. 



We represent networks as undirected and unweighted 
graphs G — {V, E) where V is the set of vertices, and E is 
the set of unweighted edges (unordered pairs of vertices). 
Multiple edges between the same pair of vertices are not 
allowed. 



The Barabasi-Albert model of scale-free 
networks 



The standard model for evolving networks with an 
emerging power-law degree distribution is the Barabasi- 
Albert model. In this model, starting from mo vertices 
and no edges, one vertex with m edges is attached it- 
eratively. The crucial ingredient is a biased selection of 
what vertex to attach to, the so called "preferential at- 
tachment:" In the process of adding edges, the probabil- 
ity Pu for a new vertex v to be attached to u is given 
by 



(1) 



*Electronie address: lioIme@tp.umu.se 



where is the degree of the vertex u. To understand the 
effect of preferential attachment, we will also investigate 
networks grown with a unbiased random attachment of 
vertices. Without the preferential attachment the net- 
works are known to have an exponential tail of the de- 
gree distribution | [To| . The time t is measured as the 
total number of added edges, which is different by factor 
m from Refs. |lOj where t is defined as the number of 
added vertices. 

It should be noted that in very large communication 
networks, such as the Internet, the users can process in- 
formation about only a subset of the whole network. How 
this affects the dynamics of network formation is inves- 
tigated in Ref. ||l^ . In the present work we neglect such 
effects and assume linear preferential attachment. 
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B. Load and capacity 

To assess the load on the vertices of a communication 
network, or any network where contact between two ver- 
tices are estabhshed through a path in the network, a 
common choice is the betweenness centrahty which 
often is seen as a vertex quantity but has a natural ex- 
tension to edges e G -E Q: 



E 



(2) 
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where (T„„, (e) is the number of geodesies between v and w 
that contains e, and cr„u, is the total number of geodesies 
between v and w. Csie.) is thus the number of geodesies 
between pairs of vertices passing e; if more than one 
geodesies exists between v and w the fraction of vertices 
containing e contributes to e's betweenness. 

In Ref. Q (see also Ref. the use of betweenness 
centrality as a load measure is given thorough motiva- 
tions. These arguments are readily generalized to the 
case of edges sensitive to overloading: Suppose that A 
is the set of pairs of vertices with established commu- 
nications through shortest paths at a given instant 
Then let A(e) denote the load oi e G E defined as the 
number of geodesies that contains e. Then we assume 
the effective load to be the average 



(3) 



where is an ensemble of A. To proceed, we restrict 
according to: 



{A : |A| = AN{N -I)}, 



(4) 



where A is constant with respect to N. This is to be 
interpreted that an element of is a set of AN{N — 1) 
pairs of distinct vertices chosen uniformly at random, and 
thus corresponds to the case where the number of estab- 
lished communication routes ending at a specific vertex 
in average increases with TV. This case can for example 
be expected in the early days of the Internet where the 
launches of new sites made the users browse a larger aver- 
age number of sites. The case where the users at average 
connects to an A'^-independent number of others is dis- 
cussed in Appendix ^ The largest approximation, when 
using the betweenness as a load measure, is probably 
that routing protocols of e.g. the Internet has implicitly 
implemented load balancing . 

To introduce overloading to the dynamics we assign 
a capacity, or maximum value, X™'^^{e) to the load, the 
same for each edge, and say that the edge e is overloaded 
if A™^'^(e) < (A(e))n. From the definition of we can 
see that our situation corresponds to having a maximum 
capacity on the betweenness centrality of the edges so 
that an edge is overloaded if Csie) > C]^'*"" (where C'^'^'^ 
is constant). If an edge is overloaded it is simply removed 



from the graph, and the betweenness recalculated, if then 
another edge becomes overloaded it is removed, and so 
on. If more than one edge is overloaded at a time we 
choose the one to remove randomly. Multiple breakdowns 
during one time step defines a "breakdown avalanche." 



C. Quantities for measuring network functionality 

To measure the network functionality we consider 
three quantities — the number of edges L, inverse geodesic 
length and the size of the largest connected subgraph 
S: For the original BA model the number of edges in- 
creases linearly as L{t) — t (i.e. one edge is added in unit 
time). But if an overload breakdown occurs in the sys- 
tem L decreases, making it a suitable simplest-possible- 
measure of the network functionality. In a functional net- 
work a large portion of the vertices should have the possi- 
bility to connect to each other. In percolation and attack 
vulnerability studies of random networks one often uses 
S to define the system as 'percolated' (or functioning), 
when the size of the largest connected subgraph S scales 
as A^ One of the characteristic features of the 

BA model networks, as well as many real-world commu- 
nication networks, is a less than algebraically increasing 
average geodesic length I. As the average geodesic length 
is infinite when the network is disconnected (as could be 
the case when an overload breakdown has occurred) we 
study the average inverse geodesic length P3[: 
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which has a finite value even for the disconnected graph 
if one defines l/ci(z;,w) = in the case that no path 
connects v and w. To monitor the fragmentation of the 
network we will also measure the number of connected 
subgraphs n. 



III. SIMULATION RESULTS 

For relative small m, typical runs are exemplified in 
Fig. |l]. For both random and preferential attachment (5) 
reaches a critical time where after the network starts to 
break down, eventually (S') reaches a steady state value. 
The breakdown develops differently in the two cases: For 
the random attachment the breakdown is relatively slow 
and the steady state value is high compared to the prefer- 
ential attachment case where large successive avalanches 
fragments the network. The other two quantities reflects 
the same behavior: While the initial vertices gets joined 
into the network increases to an early maximum. Af- 
ter the decrease corresponding to the increase of 
decreases rapidly when the network becomes fragmented. 
L shows the jagged shape, as expected correlated with 



3 



300 



7 0.15 



— 1 1 1 — 1 1 1 1 1 


: 1 1 1 1 1 1 1 1 r-p- 








r ' 




/ 1 ^ / 


. / 


/ 1 ^v^ ^ 

' l' 


/ RA 




7 PA (a) 


' (b) 


J 1000 2000 3000 4000 


1000 2000 3000 4000 


f 


f 


1 — 1 1 1 1 1 1 1 1 1 

\ 

l\ 
\\ 

. \ X 


1 1 1 l—pi 1 1 1 1 

i I u u m . 




1 1 il /"ii m ' 






'\ (c) . 


1 (d) ■ 















1.0 -J 
0.6 



1000 2000 3000 4000 1000 2000 3000 4000 
f f 

FIG. 1: The time evolution of S (a), L (b), l^^ (c), and n 
(d) for a typical run with Cs'"' = 500 and mo/2 = m = 4. 
Dashed lines represent the network grown with preferential 
attachment (PA), solid grey lines denotes curves for the runs 
with an unbiased random attachment (RA). 



that of {S}. As seen in Fig. ^(a) and (b), the disconti- 
nuity in L (in the preferential attachment case), is less 
pronounced than that in (S), so a small number of over- 
loaded edges can be enough to cause large decrease in 
(S). The reason for this behavior is that bridges (sin- 
gle edges interconnecting connected subgraphs) have a 
high betweenness and thus are prone to overloading. The 
number of connected subgraphs behaves qualitatively the 
same for random and preferential attachment. For other 
runs of the algorithm the breakdown can qualitatively 
be described as above. The averaged quantities varies 
relatively little, for example the peak-time for (S) has a 
standard deviation of ^ 3%. 

The corresponding overload case for vertices studied 
in Ref. |^ shows a similar time development with an pe- 
riod of incipient scale-freeness, an intermediate regime of 
breakdown and recovery (although the period of recovery 
is not as large for edges as for vertices), and a final break- 
down to a large-t state of disconnected clusters. One 
major difference between overload breakdown for vertices 
and edges is that the difference between random and pref- 
erential attachment is larger for edge overloading — edge 
robustness benefits more than vertex robustness from the 
geometry arising from random attachment. 

Next we investigate the m-dependence. As seen in 
Fig. 1^ the system becomes more and more robust when m 
increases. This is of course expected since with a higher 
average degree more edges shares the load, so the maxi- 
mal load can be expected to decrease. For high enough 
TO there are no avalanches, the largest connected com- 
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FIG. 2: The m dependence of the time development of 
(5') — the average size of the largest connected subgraph with 
^max _ gQQ^ g^^j _ 2y^^ ^g^-j ]-a,ndom attachment and 
(b) preferential attachment. 



ponent remains of the same size S ~ C^^^ + 1. When 
S — C'^'^^ + 1 the next edge attaching a new vertex will 
have Cb (e) = C'^'^^ + 1 , and thus be overloaded. In most 
cases this will lead to removal of the newly added edge — 
otherwise another edge has to be overloaded at the same 
time, which is decreasingly likely with increasing m. In 
Fig. |^(b) we can see one exception to this interpretation 
at TO = 6: Here (S) reaches C^^^ + 1 but starts to de- 
cay slowly at around t — 8000. As mentioned, the largest 
connected subgraph is expected to become more stable as 
TO increases. If there is an m above which (S) = C^^^ + 1 
for arbitrary large t above some to is an open question. 
Comparing Fig. ||(a) and Fig. ||(b) shows that random 
attachment and preferential attachment have similar to- 
dependence behavior — the major difference being that 
preferential attachment has a much sharper increase of 
(S) ; to be more precise the to values that does not reach 
S — C^'^^ + 1 for any t, have a lower value in the large-t 
limit. 

To get another angle of the mechanisms of the break- 
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FIG. 3: Histograms (averaged over 10'' runs) of degree. The 
parameter values are m — mo/2 — 4 and Cg^" = 500. 
(a) shows histograms for random attachment, (b) shows his- 
tograms for random attachment. The grey line in (b) is the 
function 38 k^^ illustrating the emerging power-law degree 
distribution at early times. 



downs for small m, we consider histograms of degree 
ky and betweenness CB{e). Figs. ^ and ^ shows these 
histograms both before and after the large drop in (S) 
for m = mo/2 = 4 and Cg"''' = 500. (In the ran- 
dom attachment case this drop occurs at tdrop ~ 1600, 
the corresponding value for preferential attachment is 
idrop ~ 2000.) For random attachment the difference 
between the histograms before and after the (S')-drop 
is distinctively smaller than for preferential attachment, 
just as expected from Fig. ^. The random attachment 
curves in Fig. ||(a) has a degree distribution of truncated 
exponential form both at the earlier and later times. In 
Fig. ^(a) it is exponential over two decades of P(ky), but 
falls off faster than exponentially for higher ky . Forpref- 
erential attachment the degree distributions (Fig. 0(b)) 
have a distinct difference — at t < idrop there is an emer- 
gent power-law shape of the P(/c^)-curve, whereas at 
t > idrop the shape is exponential, ~ exp(— 0.62 /c^), 
over five decades. To summarize, the degree distribu- 
tions before and after the (S')-peak illustrates the same 
behavior as the time evolution of (S) — the breakdown in 
the preferential attachment case is both faster and more 
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FIG. 4: Histograms (averaged over 10* runs) of edge between- 
ness centrality. The parameter values are — as in Fig, |^ — 
m = mo/2 — 4 and Cg^^ = 500. (a) shows histograms for 
random attachment, (b) shows histograms for random attach- 
ment. 



restructuring than in the random attachment case. 

The betweenness distributions of Fig. ^ shows a peak 
that moves to higher Cb, as t grows, until it reaches its 
maximal value at the time of the drop in (S) and starts to 
decrease. For random attachment (Fig. H (a)) the shape 
of the distribution looks qualitatively the same before 
and after the drop, but for preferential attachment (Fig. ^ 
(b)) P(Cb) ~ for betweenness smaller than the peak. 
The vertex betweenness distribution of the BA model is 
known to be strictly decreasing [T7| , which would imply 
that the low-Cs tails in Fig. ^ (b) (and most likely in 
Fig. ^ (a) as well) comes from a spread of the size of 
the largest cluster, rather than from a tail in the largest 
cluster's betweenness distribution. Another feature of 
the betweenness histograms of Fig. ^ is the smaller peaks 
at low Cb for t < idrop- These peaks corresponds to 
a sharp peak of the cluster size distribution just after 
the (S')-peak (see Fig. ||). Such smaller clusters have 
small average degree with many k^ — 1 vertices, which all 
contributes to a peak at s of the betweenness histograms. 
This explains the peak at Cb ~ 45 in the t = 5000 curve 
of Fig. I (b). 

The distribution of cluster sizes displayed in Fig. |^ 
gives some further insights: For t > idrop of the ran- 
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FIG. 5: The distribution of sizes of connected clusters at dif- 
ferent times for m — mo/2 = 4, Cg^^ = 500 and (a) random 
attachment, (b) preferential attachment. All curves repre- 
sents averages over 10* runs. To overcome noise the t > 2000 
histograms are binned for s > 75, each point being an average 
over a width of 10. 



dom attachment curves shows a bimodal distribution as 
P(s) is zero in the interval 60 < s < 290. The pref- 
erential attachment curves, in contrast, has a long tail. 
Both the large s peak for random attachment and the 
tail of preferential attachment corresponds to one single 
cluster. This is in striking contrast to the vertex over- 
load case 1^ where the network looses the unique largest 
component after the breakdown avalanches. As t evolves 
well beyond tdrop the largest component peak decreases, 
and does thus not represent a giant component (a largest 
cluster proportional to N) . The picture for both random 
and preferential attachment is thus that the system does 
not loose its unique largest cluster in a single breakdown 
avalanche — an avalanche rather results in a few isolated 
vertices or smaller clusters getting disconnected from the 
largest connected component. 

The overall picture of the time evolution of (S), L and 
(Fig.|l]), the m-dependence (Fig.||), as well as the his- 
tograms of Figs. ||, ^ and H is that for small m, avalanch- 
ing breakdowns fragments the network to a state from 
which it never recovers. For preferential attachment the 
newly fragmented network contains a single largest clus- 
ter with a well defined size, and the emergent scale-free 



degree distribution before tdrop is replaced by an expo- 
nential distribution. The breakdown for the random at- 
tachment case turns out to be less violent, and does not 
cause any major structural change. Furthermore, the dif- 
ference between the random and preferential attachment 
cases is larger for edge breakdown than for the corre- 
sponding vertex breakdown model studied in Ref. |p[. 



IV. SUMMARY AND CONCLUSIONS 

We have studied networks grown by the Barabasi- 
Albert model for networks with emergent scale-freeness 
and edges sensitive to overloading. Except the prefer- 
ential attachment defining the BA model, we also study 
an unbiased random attachment. We focus on the case 
where the number of established connections to random 
other vertices of the network scales linearly with the num- 
ber of vertices in the network. 

We find that for intermediate values of m (the number 
of edges added per vertex) the network grows like the BA 
model up to a point where is starts to break down. After 
a number of avalanching breakdowns the network reaches 
a state characterized by many disconnected clusters from 
which a giant component never re-emerges (although, in 
the preferential attachment case, there will always be one 
single largest cluster much larger than any other). If 
the growth is by random attachment, the breakdown is 
less violent with smaller avalanches and no pronounced 
structural change. For large m the steady state at large 
times is characterized by a constant largest cluster size. 

In context of real world communication networks one 
can conclude that these would benefit from being grown 
by random rather than preferential attachment (and this 
difi^erence being larger for edge overload than for vertex 
overload studied in Ref. jsj ) . In the vertex overload case 
avalanches proceeds until the network is fragmented into 
small clusters; in the edge overload problem there is still 
one large component after the breakdowns, thus we infer 
that for real-world communication networks, vertex over- 
loading is a greater threat than edge overloading, and 
congestion control in telecommunication networks 
and Internet routing protocols |^ should focus on bal- 
ancing the vertex rather than edge load. Only if the ca- 
pacity of vertices (servers etc.) grows significantly faster 
than the capacity for edges (cables etc.), edge overload 
breakdown is a potential threat for avalanching break- 
downs that is triggered by the change of load in a growing 
network. 
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FIG. 6: The time evolution of S (a), L (b), (c), and n (d) 
for a typical run in the intrinsic communication activity case. 
The model parameters are; c^^'' = 0.1 and mo/2 — m — 5. 
Dashed lines represent the network grown with preferential 
attachment (PA), solid grey lines denotes curves for the runs 
with an unbiased random attachment (RA). 



APPENDIX A: INTRINSIC COMMUNICATION 
ACTIVITY 

This paper deals mainly with the case where the aver- 
age user of a growing communication network communi- 
cates with a number of others that increases linearly with 
N. One can also imagine a case where, even though the 
network grows, the users in average communicates with a 
network size independent number of others; which is the 
topic of the present Appendix. (In Ref. ||^ this scenario 
was termed "intrinsic communication activity.") The be- 
havior of real communication networks lies, presumably, 
between these two extremes. 



1. Definitions 

To implement the situation of intrinsic communication 
activity, we modify Eq. || to: 



n' = {A: |A| = A'N} 



(Al) 



where A' is constant with respect to TV. I.e. the users has 
the iV-independent average number A' of established con- 
tacts through shortest routes. Averaging the load over 
fi' according to Eq. || gives: 



Aen' (w,w')eA 



'(e) 



A^ 

'n 



From this we see that having a constant capacity for the 
load A(e) corresponds to having a limit on CB(e) that 
increases with N. Thus we view e as overloaded if Csie) 
(where c^^'^ is constant). 



exceeds C^'"'' 



2. Results 

In the vertex overload breakdown problem, the case 
of intrinsic communication activity has a more complex 
dynamics than the extrinsic communication activity case 
(studied in the main part of the text) , with giant compo- 
nents forming only occasionally for some sets of parame- 
ter values [H . For edge overload breakdown, on the other 
hand, the dynamics of a system with intrinsic communi- 
cation activity seems very simple with no avalanching 
breakdowns and no qualitative difference between pref- 
erential and random attachment, see Fig. § We can also 
notice that the measured quantities has a power-law de- 
pendence of t. (Fig. ^ is constructed from one run with 
random and preferential attachment respectively.) For 
large times (1000 < t < 5000) the exponent a for the time 
development of the respective quantity is (in the large t 
limit): a;-i w —0.6, = 1.0, and as = an ~ 0.50 
for both (a) and (b). Initially a/-i is closer to zero, for 
^00 <t < 1000 we have ai-i « -0.5. To illustrate the 
consistency of the exponents we note that 
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(Ics) , 
(A3) 



where (^cs) is the average geodesic length for a connected 
subgraph, and d{v, w) = ii v and w are disconnected. 
This yields 



(CBie)) 



- - 1 1 ^cs 



< maxCB(e) 

vev 



(A4) 



CBie) 



(A2) 



If one assumes that {Cb{&)) c>c max^gy Cb{&) and N ^ 
n we have {Ics) oc n. Making the crude approximation 
« (^cs)~^ gives ai-i « — which holds well for 
small t. As t increases the spread in shape of the con- 
nected subgraphs becomes larger so the « (^cs)~^ 
approximation becomes worse which is seen as a slight 
increase of the slope ai-i. That the approximation 
ai~i ~ —an is rather good throughout the range of t is 
also reflected in that the average size of connected compo- 
nents N/n is never very far from S: At t = 5000 we have 
(see Fig. ^(a)) N/n « 50 and S = 57 for random attach- 
ment, and N/n « 53 and = 56 for preferential attach- 
ment. In this approximation we see that I (x n (x iV^/^ so 
the small average geodesic length is lost within the con- 
nected subgraphs. If Cg'^^ is chosen larger, so the network 
initially grows without edges breaking, there are no large 
avalanches but a crossover to the behavior seen in Fig. ra. 
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